home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / CGECO.z / CGECO
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCGGGGEEEECCCCOOOO((((3333FFFF))))                                                            CCCCGGGGEEEECCCCOOOO((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CGECO   - CGECO factors a complex matrix by Gaussian elimination and
  10.      estimates the condition of the matrix.
  11.  
  12.      If  RCOND  is not needed, CGEFA is slightly faster.  To solve  A*X = B ,
  13.      follow CGECO By CGESL.  To Compute  INVERSE(A)*C , follow CGECO by CGESL.
  14.      To compute  DETERMINANT(A) , follow CGECO by CGEDI.  To compute
  15.      INVERSE(A) , follow CGECO by CGEDI.
  16.  
  17.  
  18. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  19.       SUBROUTINE CGECO(A,LDA,N,IPVT,RCOND,Z)
  20.  
  21. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  22.      On Entry
  23.  
  24.      AAAA COMPLEX(LDA, N)
  25.         the matrix to be factored.
  26.  
  27.      LLLLDDDDAAAA INTEGER
  28.         the leading dimension of the array  A .
  29.  
  30.      NNNN INTEGER
  31.         the order of the matrix  A .  On Return
  32.  
  33.      AAAA an upper triangular matrix and the multipliers
  34.         which were used to obtain it.
  35.         The factorization can be written  A = L*U  where
  36.         L  is a product of permutation and unit lower
  37.         triangular matrices and  U  is upper triangular.
  38.  
  39.      IIIIPPPPVVVVTTTT INTEGER(N)
  40.         an integer vector of pivot indices.
  41.  
  42.      RRRRCCCCOOOONNNNDDDD REAL
  43.         an estimate of the reciprocal condition of  A .
  44.         For the system  A*X = B , relative perturbations
  45.         in  A  and  B  of size  EPSILON  may cause
  46.         relative perturbations in  X  of size  EPSILON/RCOND .
  47.         If  RCOND  is so small that the logical expression
  48.         1.0 + RCOND .EQ. 1.0
  49.         is true, then  A  may be singular to working
  50.         precision.  In particular,  RCOND  is zero  if
  51.         exact singularity is detected or the estimate
  52.         underflows.
  53.  
  54.      ZZZZ COMPLEX(N)
  55.         a work vector whose contents are usually unimportant.
  56.         If  A  is close to a singular matrix, then  Z  is
  57.         an approximate null vector in the sense that
  58.         NORM(A*Z) = RCOND*NORM(A)*NORM(Z) .  LINPACK.  This version dated
  59.      08/14/78 .  Cleve Moler, University of New Mexico, Argonne National Lab.
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCGGGGEEEECCCCOOOO((((3333FFFF))))                                                            CCCCGGGGEEEECCCCOOOO((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      Subroutines and Functions LINPACK CGEFA BLAS CAXPY,CDOTC,CSSCAL,SCASUM
  75.      Fortran ABS,AIMAG,AMAX1,CMPLX,CONJG,REAL
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.